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Weakly nonlinear hexagon convection patterns coupled to mean flow are investi- 
gated within the framework of coupled Ginzburg-Landau equations. The equations 
are in particular relevant for non-Boussinesq Rayleigh-Benard convection at low 
Prandtl numbers. The mean flow is found to (1) affect only one of the two long- 
wave phase modes of the hexagons and (2) suppress the mixing between the two 
phase modes. As a consequence, for small Prandtl numbers the transverse and the 
Qh| longitudinal phase instability occur in sufficiently distinct parameter regimes that 
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Abstract 



they can be studied separately. Through the formation of penta-hepta defects, they 

lead to different types of transient disordered states. The results for the dynamics 

of the penta-hepta defects shed light on the persistence of grain boundaries in such 

jy-j ■ disordered states. 
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1 Introduction 
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Roll patterns in Rayleigh-Benard convection of a fluid layer heated from below 
have been explored extensively over the years as a paradigmatic system to 
study the succession of transitions from ordered to disordered and eventually 
turbulent states (for a recent review see [1]). Substantial theoretical effort has 
been devoted to identifying all the linear instabilities that straight roll patterns 
can undergo (e.g. [2]). Depending on the wavenumber of the roll pattern, 
the temperature difference across the layer (Rayleigh number), and the fluid 
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properties (Prandtl number) a host of instabilities has been identified. The 
main instabilities being the Eckhaus, zig-zag, cross-roll, oscillatory, and skew- 
varicose instabilities. They limit the band of stable wavenumbers of straight 
rolls (Busse balloon). It turns out that the oscillatory and the skew- varicose 
instabilities depend sensitively on the Prandtl number and are relevant only 
if the Prandtl number is sufficiently small [2] . 

For small convection amplitudes a weakly nonlinear description in terms of a 
Newell- Whitehead-Segel equation [3,4] would be expected to be sufficient at 
least for almost straight roll patterns. However, it is found that this is true 
only in the limit of large Prandtl numbers. It has been shown that slightly 
curved roll patterns drive a mean flow. Because of the incompressibility of 
the fluid it induces a non-local coupling of the rolls, which is not contained 
in the Newell- Whitehead-Segel equation [5]. The strength of the mean flow 
increases with decreasing Prandtl number. For small Prandtl numbers the 
Newell- Whitehead-Segel equation has therefore been extended to include an 
equation for the vertical vorticity characterizing the mean flow [5]. 

Experimentally, the mean flow has been observed in experiments in a rela- 
tively small circular container [6]. Since the usual boundary conditions re- 
quire the convection rolls to be essentially perpendicular to the side-walls the 
rolls become strongly curved in this geometry. The resulting mean flow has 
been identified as the driving force for the observed persistent creation and 
annihilation of dislocations [7-9] . 

Arguably the most interesting state that is due to the mean flow associated 
with small Prandtl numbers is the spiral-defect chaos observed in large-aspect 
ratio experiments on thin gas layers [10]. It is characterized by the appearance 
of various types of defects in the pattern with small rotating spirals being the 
ones that are visually most striking. Spiral-defect chaos does not arise from a 
linear instability of the straight roll pattern and, in fact, bistability of straight 
roll patterns and spiral-defect chaos has been observed experimentally [11]. 
The onset of spiral-defect chaos depends strongly on the Prandtl number, 
indicating that mean flows play an essential role in maintaining this state 
[12-15]. 

Motivated by the strong impact of mean flows on roll convection patterns we 
consider in this paper the effect of such flows on the stability and dynamics 
of hexagonal patterns. Hexagonal patterns are commonly found in spatially 
extended non-equilibrium systems such as non-Boussinesq Rayleigh-Benard 
convection (e.g. [16]), Marangoni convection (e.g. [17]), Turing structures in 
chemical systems [18], crystal growth (e.g. [19]), and surface waves on ver- 
tically vibrated liquid or granular layers (Faraday experiment) [20,21]. Not 
in all of these systems mean flows of the type discussed above arise. Clearly, 
they are relevant for the small Prandtl numbers in gas convection in very 



thin layers [10]. Interestingly, the skewed- varicose instability and a (transient) 
state similar to spiral-defect chaos have been observed also in vertically vi- 
brated granular layers [22,23]. Since in convection at low Prandtl numbers 
they are a signature of the importance of mean flow, mean flow may also be 
relevant in vertically vibrated fluids. We focus in this paper on patterns aris- 
ing from a steady bifurcation, as is the case in convection. Immediately above 
onset, parametrically excited standing waves like those arising in the Faraday 
experiment behave in many respects like patterns that are due to a steady 
bifurcation (e.g. [24]). Our results may therefore also be relevant for hexagon 
patterns in suitable Faraday experiments. 

In the absence of mean flows the stability of hexagonal patterns has been 
studied in detail in the weakly nonlinear regime. Starting from three coupled 
Ginzburg-Landau-type equations for the amplitudes of the rolls that make up 
the hexagonal pattern two coupled phase equations have been derived that 
describe the dynamics of long-wave deformations of the hexagon patterns [25- 
28]. In these theoretical analyses two types of long- wave instabilities have been 
identified, a longitudinal and a transverse mode, with the longitudinal mode 
being relevant only for Rayleigh numbers very close to the saddle-node bifur- 
cation at which the hexagons first appear. In a more detailed analysis also 
perturbations with arbitrary wavenumber and simulations of the nonlinear 
evolution ensuing from the instabilities have been included [29]. The insta- 
bilities typically lead to the formation of penta-hepta defects [30,31]. Their 
dynamics and interaction has also been studied in the absence of mean flows 
[32,33]. In the strongly nonlinear regime long- wave perturbations are still cap- 
tured by phase equations [28]. Specifically for convection at large Prandtl num- 
bers driven by a combination of buoyancy and surface tension the stability of 
hexagons and their dynamics has been investigated in [34]. Experimentally, 
the stability of hexagonal convection patterns has not been studied in detail. 
A contributing factor has been that controlled changes in the wavenumber 
of the pattern are considerably more difficult to effect than in systems like 
Taylor-vortex flow, where detailed agreement between experiment and theory 
has been achieved (e.g. [35,36]). Recently, however, it has been possible to 
use localized heating as a printing technique for convection patterns [37,38]. 
This will make detailed analyses of the wavenumber- dependent instabilities of 
hexagon patterns experimentally accessible. 

The goal of this study is to investigate the linear stability of hexagons coupled 
to a two-dimensional mean flow driven by modulations of the pattern and 
the nonlinear evolution arising from the instabilities. An example of the flow 
generated by two typical patterns is shown in fig.la,b. In fig. la a periodic 
hexagon pattern is sinusoidally modulated on a length scale five times the 
hexagon size and fig. lb shows the mean flow generated by a penta-hepta defect. 
The mean flow is found to couple only to one of the two long-wave modes. As 
a consequence, for sufficiently small Prandtl numbers each of the two long- 
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Fig. 1. Illustration of the mean flow due to modulation of the hexagon envelope 
(panel(a)) or due to a penta-hepta defect (panel(b)). For visual clarity the wave- 
length of the hexagon pattern has been chosen larger (relative to the modulation 
wavelength) than is appropriate for the weakly nonlinear theory presented in this 
paper. 

wave instabilities dominates in a separate, experimentally accessible parameter 
regime. Both instabilities induce the formation of defects. We study their 
motion and briefly touch upon its impact on disordered patterns and grain 
boundaries. 

The paper is organized as follows: first we formulate the problem extending 
previous work on the mean flow generated by roll convection [39,40]. We then 
study the stability of the hexagonal pattern with respect to long-wave pertur- 
bations (sec. 3) and with respect to general side-band perturbations (sec. 4). In 
section 5 the nonlinear evolution of the instabilities including the formation, 
dynamics, and annihilation of defects is investigated. There, we also study the 
transition from hexagonal to roll patterns triggered by the side-band instabil- 
ities. 



2 Amplitude equations 



The interaction between weakly nonlinear convection rolls and the mean flow 
generated by them has been investigated for Rayleigh-Benard convection with 
both no-slip and stress-free boundary conditions at the top and bottom plate 
[39,40]. The relevant amplitude equations were derived and used to determine 
the stability of rolls with respect to side-band instabilities. Expanding the 
fluid velocity in terms of the complex amplitude A of the convection rolls and 
the real streamfunction Q for the two-dimensional mean flow as (e <C 1), 



v(x, y, z, i) = eA(x, y, t)e^v q (z) + e 2 { V x (-Qz)f(z)} + h.o.t, (1) 
the extended Ginzburg-Landau equation for A and Q is given by 



d t A = (fj, + (d x - i\d 2 y f - \A\ 2 )A - i Sl Ad y Q + h.o.t. 
M(Q)=2 qi (d x - i\d 2 y )d y \A\ 2 + h.o.t. 



(2) 
(3) 



where A = e/2q c with e the supercriticality parameter and q c the critical 
wavenumber. Here M(Q) = V 2 Q for no-slip boundary conditions [39] and 
M(Q) = (?2V 2 9 t — V 4 )Q for the stress-free case [40]. It is worth noting that 
in the no-slip case the mean flow is not an additional dynamical variable and 
could in principle be adiabatically eliminated which would lead to a non- 
local equation for A. The slow time t is given by t — e 2 t. The non-isotropic 
Newell- Whitehead-Segel operator (d x — i\d 2 ) 2 reflects the anisotropic scaling 
of the slow spatial scales x = ex and y = t l ^ 2 y. To leading order, the mean 
flow is driven by variations in the magnitude of the convective amplitude and 
couples back through the advective term is\Ad y Q = is\AV x , where V x is the 
re-component of the mean flow. At higher orders a term of the form V ■ VA 
would arise as well. 

For the description of hexagonal convection patterns coupled to the mean 
flow the treatment has to be extended to include rolls in other directions 
than the rr-axis. Using a coordinate transformation, it is straightforward to 
generalize the coupling term Ad y Q for rolls in the x-direction n\ = (1,0) 
to rolls in the other two hexagonal directions fi2 = (— 1/2, y / 3/2) and n 3 = 
(—1/2, — y / 3/2). The amplitudes A^ % = 1,2,3, correspond then to rolls with 
wave vector q c fii. For the generalization of the advection term iAd y Q in (2) it 
is convenient to introduce the unit vector fj perpendicular to h^. To complete 
the extension of eqs.(2,3) from roll to hexagonal patterns we also add non- 
Boussinesq quadratic terms. For no-slip boundary conditions, the minimal set 
of equations describing hexagons coupled to a mean flow then reads 



dtAj 



/i+(n j -V) 2 -|A i | 2 -z/(|A j _i| 
A)_ l A) +l ~iPA ] {f r V)Q, 
V 2 Q = -{-2d xy \A 1 \ 2 + 
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(4) 



(5) 



with cyclic permutation on j. It is not clear how the different scaling of the 
longitudinal spatial variable in the direction of hi and of the transverse variable 
in the direction of fi can be implemented in the hexagonal case where these 
directions are different for each of the three hexagon modes. We therefore 



use an isotropic scaling, x = ex and y = ey, which renders the derivatives 
i\(fj ■ V) 2 higher order in e. This can lead to degeneracies in certain growth 
rates like that of the zig-zag instability (e.g. [41,42]). Then higher-order terms 
have to be re-introduced. For the hexagonal patterns this is, however, not the 
case. Note that while there is no direct damping of transverse variations, each 
of the Aj is still allowed to vary in its respective transverse direction. This 
variation induces a variation in the longitudinal direction of the other two 
modes, in particular through the resonant interaction term A*__ X A* +1) which 
then induces damping. 

The cubic coefficient v decreases monotonically with Prandtl number from 
v ~ 2.0 for Pr = 0.5 to v ~ 1.5 for Pr ~ 10 [43,16]. For the numerical results 
in this paper we fix v = 2. Note that the mean flow amplitude Q is of higher 
order than the hexagonal amplitudes, in fact, Q oc 0(A 2 ). Nevertheless, the 
coupling term iA(fj ■ V)Q is retained because \f3\ can be much larger than 
one for low Prandtl number. For example, f3 ~ —18 for Prandtl number ~ 0.1 
[39]. Furthermore, the mean flow introduces qualitatively new effects since it 
provides a non-local coupling between the convection amplitudes that is due 
to the incompressibility of the fluid. 

In the absence of the mean flow, there exists a Lyapunov functional T for (4), 
which guarantees that only stationary patterns are possible as t — > oo. Simple 
stationary solutions are: 



(i) rolls A 3 = {ii-q 2 fl\ A j±1 = 0, (6) 



A A A l±V /l + 4 ^-^)( 1 + 2Z/ ) 

v n) hexagons A\ = A 2 = A 3 = , (7) 



2(1 + 2i/) 



'iii) mixed modes A~ = , A~ 



I pi- q 2 - A\ 



v-Y J ^ V 1 + v 

As usual, the hexagon solution given by equation (7) with a minus sign in 
front of the square root is not stable. Thus, in the following we focus on the 
hexagon solution with the plus sign. 

Since the mean flow couples to hexagons only through a modulation in the 
amplitudes, the stationary states (6-8) with constant amplitudes are also so- 
lutions to equations (4,5). However, their stability is modified by the coupling 
to the mean flow. In the case of rolls it can lead to an oscillatory instabil- 
ity (oscillatory skew- varicose instability), implying that there is no Lyapunov 
functional when taking the mean flow into account. 



3 Linear Stability Analysis: Long-wave Approximation 



We conduct a linear stability analysis of the hexagonal pattern (7) in the 
long- wave limit following the procedures in [26] . The amplitudes are perturbed 
around a hexagon pattern, 



Aj = RoeW<*>v)(l + r 3 + ify), j = 1, 2, 3, 



(9) 



where R = (l + yl + A(fj, - q 2 )(l + 2u)j /2(1 + 2v) is the amplitude of the 
stationary solution with reduced wavenumber q, and Tj and <f)j are the am- 
plitude and phase perturbations, respectively. We also define u and v as in 



u = R 2 {1 - u) + Ro, v = 2R 2 (l + 2v)-R , 



(10) 



where v — corresponds to the saddle-node bifurcation at which the hexagons 
come into existence. It is equivalent to 



Vsn = q 



1 1 



41 + 2i/ 



At u = the hexagons become unstable to the mixed mode solution 
equivalent to 



(11) 
). It is 
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Both u and v have to be greater than zero for hexagons to be stable. 



(12) 



We first rewrite (5) in terms of the perturbation amplitudes. The mean flow 
then satisfies the linearized equation 



V 2 Q = -2Rl{-2d xy n + 



V3. 



-d 2 
v 



dl) + a 



Xjl 



r 2 + 



2 y y 



d 2 ) + a 



Xjl 



rz}- (13) 



Substituting the perturbed amplitudes (9) into (4), we obtain 



d T n = -An + Ai(r 2 + r 3 ) + {hi ■ V) 2 ri - 2q{hi ■ V)0i, 
d T r 2 = -Ar 2 + Ai(r 3 + n) + (n 2 • V) 2 r 2 - 2g(n 2 • V)0 2 , 
,9 T r3 = -Ar 3 + Ai(n + r 2 ) + (n 3 • V) 2 r 3 - 2g(n 3 • V)0 3 , 



(14) 
(15) 
(16) 
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drfc = -i? $ + 2q{h 1 • V)n + (ni • V) 2 0i - SifyQ, (17) 

9t0 2 = -#0$ + 2g(n 2 • V)r 2 + (n 2 • V) 2 2 - Sl ( ^-d x - \d y ) Q, (18) 



<9 T 3 = -i? $ + 2g(n 3 • V)r 3 + (n 3 • V) 2 3 - Si \^d x - ld y ) Q, (19) 



where A = (1 + 2R )R , A x = (1 - 2uR )R . 

In this section we focus on the limit of long-wave perturbations and introduce 
super-slow scales X = 5x, Y = Sy, and T = 5 2 t and a slow gradient, V = 5V 
(5 <C 1). In this limit the amplitude perturbations Tj and the global phase 
$ = 0! + 2 + 3 are slaved to the two slow modes (j) x = — (0 2 + 03) and 
(f) y = (02 — 03)/v^3, which arise from the translation symmetry of the system. 
Adiabatic elimination reduces then (14-19) to 



n=r 3 - ^((nx • V)0x - (n 3 • V)0 3 ), (20) 

u 

r 2 = r 3 -^((n 2 -V)0 2 -(n 3 -V)0 3 ), (21) 

u 

r 3 = -^(n s • V)0 3 - g(2 ^~ ^ fork • V)0, - 3(n 3 • V)0 3 ) , (22) 

3 3 

3R ^ = J2(n l -V) 2 (j) i + 2qJ2(n i V)n, (23) 



and an evolution equation for the phase vector 



'xi fy)i 



dt$= D ± V 2 + (£>,, - £> ± )V(V • 0) + M^V x (V x 0). (24) 

2m 



The last term of the phase equation (24) arises from the coupling to the mean 
flow, which is driven by the curl of the phase, 

Q = J_3l { _ gz .vx0). (25) 

lu 



In (24) D» and D±_ are the longitudinal and transverse phase diffusion coeffi- 
cients for the hexagons in the absence of a mean flow [25,27,44]. 

Even in the presence of the mean flow, equation (24) allows a decomposition 
of the phase vector into a longitudinal (curl-free) mode <pi satisfying 



e z • V x ^ = 0, (26) 




5 10 15 

Prandtl number Pr 
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Fig. 2. The coupling strength (3 for rolls. Solid line from [39]. Dashed line is the 
dependence conjectured in [45] with the free constant chosen to match the solid line 
at Pr = 0. 

which has growth rate 
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(27) 



and a transverse (divergence-free) mode <p t satisfying 



V- 



[2%) 



with growth rate 



0t 



k* 

2 



u 



3qf3R t 



a 



(29) 



To this order the growth rates do not depend on the orientation of the per- 
turbation wave vector k = (k, I) but only on its magnitude k = \Jk 2 + P. 

As expected from equation (25), the mean flow only affects the transverse 
mode. It shifts the associated stability boundary of the hexagons and makes 
it asymmetric in q. The sign of (3 determines whether the mean flow desta- 
bilizes the transverse mode for positive or for negative wavenumbers. For roll 
convection with no-slip boundary conditions (5 is always negative [39] and its 
dependence on the Prandtl number is shown in figure 2. We note that /3 ~ — 1 
for gases (unity Prandtl numbers), and (3 ~ —0.1 for water (Prandtl number 
~ 10). 
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Fig. 3. Summary of the stability boundaries for infinite Prandtl numbers {(3 = 0) 
showing the saddle- node bifurcation (thin solid line), transition of hexagons to 
mixed-mode solution (thin dashed), longitudinal mode (thick solid), transverse 
mode (thick dotted), equal-energy line for hexagons and rolls (thick dashed), 
cross-over from longitudinal to transverse mode (thick dashed-dotted). 

The stability boundaries for infinite Prandtl number (no mean flow) are re- 
viewed in figure 3. Hexagons exist above the thin solid line (v — 0). The thick 
solid line and the dotted line are the stability boundaries for the longitudi- 
nal and transverse modes, respectively. Along the thick dashed line, rolls and 
hexagons have the same value of the Lyapunov functional. Hexagons become 
unstable to rolls (via the unstable mixed mode (8)) at the thin dashed line 
(u — 0). Note that since we have scaled the coefficient of the quadratic nonlin- 
earity in (4,5) to 1, this transition occurs at /i = 4. For weak non-Boussinesq 
effects this corresponds to very small convection amplitudes. As discussed be- 
low (see also sec. 4), the decomposition into transverse and longitudinal mode 
does not hold beyond the long- wave limit. It is replaced by a separation into a 
'wide-splitting' and a 'narrow-splitting' mode [29]. The dashed-dotted line sep- 
arates the regions in which the wide-splitting and the narrow-splitting mode 
represents the mode with maximal growth rate, respectively. In view of the 
effect of the mean flow it is worth noting that the regime in which the longi- 
tudinal mode is the relevant destabilizing mode is typically extremely small 
(below dashed-dotted line in fig.3) and it is very difficult to investigate that 
instability experimentally. In fact, even in numerical simulations of (4) (with 
(5 = 0) we found it difficult to separate the dynamics of the two modes (see 
also [29]). Note that due to the scaling of the quadratic coefficient in (4,5) the 
usually small range over which hexagons can be observed extends from \x rj 
to // ~ 1.5. 

Figure 4a,b shows how the stability boundaries are altered by the mean flow 
for (5 = —0.2 and (5 = —3.0, respectively. Figure 4(c) shows an enlargement 
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Fig. 4. Stability boundaries of hexagons coupled to a two-dimensional mean flow. 
The thick solid line is the stability limit given by the longitudinal mode. The thick 
broken lines are for the transverse mode. Thin lines are as in fig. 3. Panel (a) is for 
f3 = —0.2 (Pr ~ 5) and panel (b) is for (3 = —3 (Pr ~ 0.5). Panel(c): enlargement 
of the stability boundaries for small fi showing results for f3 = 0, f3 = —0.2, and 
(3 = -3. 

of figure 4a,b for small \x. As discussed before, the mean flow affects only 
the transverse mode. The corresponding stability limits are denoted by thick 
dotted, dashed and dashed-dotted lines for f3 — 0, (5 — —0.2, and (5 = —3.0, 
respectively. The thick solid line gives the stability limit due to the longitudinal 
mode, which is independent of (3. For all values of /3 it is the relevant instability 
immediately above the saddle-node bifurcation at which the hexagons first 
arise (cf. fig.4c). For somewhat larger values of /x the transverse mode becomes 
relevant as well. In the physically relevant case, f3 < 0, it is destabilized by 
the mean flow for positive q strongly reducing the stability balloon there. 
While for large Prandtl numbers the cross-over from the longitudinal to the 
transverse mode is almost imperceptible in the stability limit itself, it occurs 
very suddenly when the mean flow is important. This behavior has also been 
found in fully nonlinear stability calculations of hexagonal convection in non- 
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Boussinesq convection [46] . For q < the transverse mode is stabilized by the 
mean flow. This allows the longitudinal mode to become important again for 
larger //. In fact, for sufficiently large /3 the transverse mode does not contribute 
to the stability limits at low q for any value of \x (cf. dashed-dotted line for 
f3 = —3 in fig.4b). Thus, for small Prandtl numbers the mean flow makes 
it possible to investigate the two instabilities separately. In sec. 5.1 below we 
discuss the different nonlinear behavior arising from the two instabilities. 

The leading-order phase equation (24) is isotropic. Only at higher order in 
the gradients this symmetry is reduced to that of the underlying hexagon 
pattern. Similarly, at higher orders (i.e. for larger magnitudes of the perturba- 
tion wavevector) the eigenmodes of the phase equation are not expected to be 
purely transverse or longitudinal any more. In sec. 4 we determine the charac- 
ter of the destabilizing modes numerically for finite perturbation wavevectors, 
but only for specific values of the parameters. To obtain more general results 
for the influence of the mean flow on the destabilizing modes we derive here 
the phase equation to fifth order. 

If only one of the two modes is unstable the coupled phase equations can 
be reduced to a single equation for the unstable mode [26]. It is obtained by 
eliminating the stable mode adiabatically. Since both modes are long-wave 
modes satisfying diffusion-like equations the elimination introduces non-local 
terms in the equation for the unstable mode. To demonstrate more explic- 
itly the breaking of the isotropy and the mixing of the two modes we derive 
the coupled fifth-order equations for both components of the phase vector <fi. 
This can be done systematically in the vicinity of the codimension-two point 
(q^ ct \ /i^) where both oi and at are zero to cubic order (cf. (27,29)). There the 
orientation of the eigenmodes depends solely on the fifth-order terms. Away 
from the codimension-two point it is determined by the competition between 
the cubic and the fifth-order terms. 

To facilitate the analysis we restrict to cases of weak mean flow and expand in 
|/3| <C 1. To cubic order, the codimension-two point is given by q^ = q + (3q 1 
and Rq = R 00 + (3R m with 



V^+2 „ _ 3(i/+l) _ 1 _ 3y^7+2 



Here we consider the magnitude Rq as the control parameter instead of fi. To 
first order in /3, the fifth-order, linear phase equations at the codimension-two 
point read as follows: 
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dtV 2 ^ = -^ryV 2 [{Adld y + 12d x d 3 y )<f> y + (ll^ 4 + 3d 4 + 6^)0, 
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Fig. 5. Panel (a): angles i/j and 9. Panel (b): Angle deviation Sip for the two eigen- 
vectors at the codimension-two point as a function of the angle 6 between k and n\ 
for v = 2. The thin solid line is for SipQ ' , and the thick dash-dotted line is 5tp\ 



and the thick dashed line is 5ip\ 



(2) 



2(u + l) 



+p w(u + i) 2 ^( 3(4 + 13z/)<9 ^ + 6(8 + 3u)d * d y + 9(4 + 3z/ )^4 5 ) h 

+ (3(11 + 7v)dl + 3(17 + 26v)d*d 2 y + 9(3 + v)d 2 x d* + 9d 6 y ) <p x ] , (31) 



d t V 2 



v 



2(i/ + 1) 



(<9 4 + 9«9 4 + 18«W)0, + (4$jfo + 12^ 3 )<A 



(3Sg + 3(19 + 13u)d 4 x d 2 y + 9(9 + 2u)d 2 x d 4 y + 21d\ 



+/? 16(z/+l) 2 

+ (3(4 + lv)d\d y + 3(8 + Yiv)dld\ + 9(4 + u)d x tf) <p x 



(32) 



A brief summary of the derivation and the full nonlinear phase equations (for 
/3 — 0) are included in the appendix. Similar to the nonlinear phase equations 
derived in [26], eqs. (31,32) are non-local. Here, however, the reason for this is 
the nonlocal effect of the mean flow: for j3 = the equations become local. 

The mixing between the transverse and the longitudinal modes can be quan- 
tified by the angle ip between the phase vector and the wave vector k of the 
perturbation eigenmode, 



ip = arctan 



,e z ■ V x 
' V ■$ 



(33) 



Figure 5 (a) illustrates the relative positions of <p and k with respect to n\. 
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(1 2) 

Denoting the phase angle for the pure modes by if>Q ■> the pure transverse 
mode has ipQ — tt/2 while if>Q = for the pure longitudinal mode. From 
eqs. (31,32) we find that one eigenvector, 1; is almost transverse (if> ~ 7r/2) 
whereas the other, <f) 2 , is almost longitudinal (if) ~ 0). For each mode we 
expand the deviation from if>Q in (3 at the codimension-two point: 

^(i,2) = ^,(1,2) _ ^,(1,2) = ^(i.2) (l/> 0) + ps$' 2 \v, 9) + 0(/3 2 ), (34) 



where is the angle between k and ni. It turns out that in the absence of mean 
flow the two eigenvectors remain orthogonal to each other to fifth-order and 
therefore SipQ = 5ipQ . The deviation 5if) ' depends very little on v. The 

(1 2) 

contributions 8if)\ ' from the mean flow decrease with increasing v without 

Q 2} 

changing sign. Fig. 5 shows as an example the contributions 5i/jq { to 5ip for 

(1 2) 

v — 2. The thin solid lines are 8if>Q ' for the two eigenvectors. The thick 
dashed-dotted line is 8if>[ for 0i and the thick dashed line is for 02- Note 
that, because 5ip\ ' and SipQ ' have the same signs and because /3 is negative, 
the mean flow always suppresses the mixing between the two modes. Thus, 
the phase perturbations are closer to the pure transverse/longitudinal modes 
for small Prandtl number. This is also found for parameters away from the 
codimension-two points (see fig. 6 below). 



4 General Linear Stability Analysis 



In this section we present results from solving the linear-stability equations 
(14 - 19) numerically for arbitrary perturbation wavenumbers. We first review 
and discuss the results for infinite Prandtl number [25,29,27]. The long- wave 
analysis of sec. 3 revealed two phase instabilities associated with the longi- 
tudinal and the transverse phase mode, respectively. For finite perturbation 
wavenumbers these two modes become mixed and the eigenmodes of the linear 
stability problem are not strictly longitudinal or transverse any more. Instead 
they can be characterized by the angle 9 between the perturbation wavevector 
of the fastest growing modes and n\. For the 'wide-splitting' modes 9 = tvk/3, 
whereas for the 'narrow-splitting' modes 9 = ir/6+nir/3 with n an integer [29]. 
The terminology refers to the fact that the angle between the relevant most 
unstable wide-splitting modes is 2-7r/3, whereas for the narrow-splitting modes 
it is only 7r/3 [29]. As the perturbation wavenumber goes to the narrow- 
splitting mode turns into the longitudinal mode and the wide-splitting mode 
into the transverse phase mode. This is illustrated in fig. 6, which shows the 
angle if) (eq.(33)) as a function of the angle 9 between k and Hi. The thin 
lines refer to the case /3 — discussed here. The deviation of the destabilizing 
modes from these orientations increases with |k| (cf. fig.5). 
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Fig. 6. Angle ip (between the phase vector cj> and the perturbation wave vector k) 
as a function of the angle 8 (between k and Hi) for \x = 0.5 and q = 0.5. The thin 
lines are for coupling strength = 0, and (3 = —0.2 for the thick lines. The solid 
lines are for k = 0, the dashed lines for k = 0.192, and the dashed-dotted lines for 
k = 0.224. 
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Fig. 7. Cross-over from narrow-splitting to wide-splitting mode. Along the 
dashed-dotted (dashed-double-dotted) line the narrow- and wide-splitting modes 
have equal growth rates for /? = (/? = —0.2). 

Numerical simulations show that the nonlinear evolution of the narrow- and of 
the wide-splitting modes leads to qualitatively different transient patterns with 
the wide-splitting mode inducing noticeably more grain boundaries between 
hexagon patterns of different orientation and consequently a larger number 
of defects than the narrow-splitting mode [29]. However, as fig. 3 already sug- 
gests, the range of parameters over which the longitudinal/narrow-splitting 
mode dominates over the transverse/wide-splitting mode is very small and it 
is very difficult to separate the two modes in numerical simulations [29]. This 
is illustrated further in fig.7 and in fig. 8. Fig. 7 gives the cross-over line from 
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Fig. 8. Results of full linear stability analysis with (5 = 0. Panel (a): maximal 
growth rate as a function of jjl for q = 0.6, 0.45 and 0.35 (curves 1 to 3). Panel (b): 
wavenumber of the fastest growing mode for the corresponding q (curves 1 to 3) . 

the narrow-splitting to the wide-splitting mode, i.e. the line along which the 
growth rates of both modes are equal. Fig.8a shows the dependence of the 
maximal growth rate of the dominating mode as a function of the control pa- 
rameter fj, for three values of the wavenumber. Fig. 8b shows the dependence 
of the wavenumber of the perturbation mode with maximal growth rate. As 
expected for a long-wave instability, the magnitude of the wavenumber of the 
fastest growing mode goes to as the stability limits are reached at fi ~ 0.6 
and /i ~ 1.2 for q = 0.35. Although the relevant perturbation wavenumber 
jumps by almost a factor of 2 across the transition from the narrow-splitting 
to the wide-splitting mode (cf. the jump at /i ~ 0.4 for q = 0.6), it turns out 
that in simulations this difference is very hard to identify because the growth 
rates of the two modes are so similar and the narrow-splitting mode dominates 
only over such a small range of ft. 

The mean flow introduces a number of modifications in the stability behavior. 
The linear operator £ becomes non-hermitian and allows complex eigenval- 
ues. It turns out, however, that in our calculations complex eigenvalues appear 
only near the stability boundary of the mixed mode (/i = Hmm) and beyond 
the long- wave stability limit given by the longitudinal mode (cf. fig. 4). In nu- 
merical simulations, the transient oscillations were always accompanied by the 
formation of defects and an eventual transition to rolls. Possibly, the mean flow 
induces more relevant oscillatory instabilities in combination with nonlinear 
gradient terms, which also constitute non-potential terms [47,48]. 

The mixing between the longitudinal and the transverse mode for finite values 
of k persists in the presence of mean flow. It is, however, strongly suppressed, 
consistent with our analytical results for the vicinity of the co dimension-two 
point (cf. eq.(34)). This is illustrated in fig. 6, where the thick lines give the 
angle ip of the phase vector for (3 = —0.2. The difference between the growth 
rates of the narrow- and of the wide-splitting modes is still quite small for 
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Fig. 9. Results of full linear stability analysis with q = ±0.6. Solid lines are for 
q = 0.6 and dashed lines are for q = —0.6. Panel (a): maximal growth rate as a 
function of \x. Panel (b): wavenumber for the fastest growing mode. From curves 
to 3 = 0, -0.1, -0.2, -0.3, respectively. 

small values of /i, as it is in the case without mean flow. This is illustrated in 
fig. 9, which shows the maximal growth rates and the associated wavenumbers 
for q = +0.6 and for q = —0.6 using solid and dashed lines, respectively, for 
different coupling strengths. For larger values of /x, however, even moderate 
values of the coupling enhance the difference noticeably. Thus, for /3q < the 
wide-splitting mode becomes much more dominant. For j3q > the narrow- 
splitting mode, which in the absence of mean flow is only relevant in a very 
narrow range of the control parameter near the saddle-node bifurcation, be- 
comes dominant again above a critical value of /i. This reflects the asymmetry 
of the stability boundaries induced by the mean flow as shown in fig.4a,b. 
The critical value of \x decreases with decreasing Prandtl number, and below 
Prandtl number ~ 2.6 (/3 = —0.4), the transverse mode is preempted by the 
longitudinal mode over the whole range of ji. This suggests again that for suf- 
ficiently small Prandtl numbers it should be experimentally possible to study 
the two instability modes separately. 



5 Numerical simulation 



We now present results from numerical simulations of eqs.(4,5) using a parallel 
pseudo-spectral code. In sec. 5.1 we compare the nonlinear evolution of the 
transverse and longitudinal side-band instabilities. We then discuss the impact 
of the mean flow on the competition between rolls and hexagons. In sec. 5. 2 we 
study the motion of penta-hepta defects in the presence of mean flow. These 
results shed some light on the persistence of grain boundaries in disordered 
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hexagons patterns. 

5.1 Nonlinear Evolution of the Side-band Instabilities 

Transverse vs. longitudinal modes. 

Without mean flow the nonlinear evolution of the two side-band instabilities of 
hexagon patterns has been studied numerically in some detail by Sushchik and 
Tsimring [29] . They find that both instabilities lead to the formation of penta- 
hepta defects and subsequently to disordered hexagon patterns. When only the 
wide-splitting perturbations are growing the emerging patterns have many 
grain boundaries separating domains of hexagons with different orientation. 
These disordered patterns persisted for a long time. When the narrow-splitting 
instability is dominant only few such domains appear and a regular hexagon 
pattern with a different wavenumber is restored within a relatively short time. 
Sushchik and Tsimring note that numerical simulations in the regime in which 
the narrow-splitting mode is the only destabilizing mode are very difficult and 
only present results for parameter values for which also the wide-splitting 
mode is destabilizing albeit to a lesser degree. These runs correspond to a 
quench deep into the unstable regime beyond the dashed-dotted line in fig. 3 
(q = 0.56 and \x = 0.6). The authors associate the larger number of grain 
boundaries in the first case with the large angle between the components of 
the wide-splitting mode. A detailed comparison is difficult since no simulations 
are available in which only the narrow-splitting mode is active. 

As discussed in sec. 3 and sec. 4, the mean flow couples only to the transverse/wide- 
splitting instability and for sufficiently small Prandtl number suppresses it 
completely for f3q > 0. This makes it possible to compare the two instabilities 
in detail under comparable conditions. Figs. 10-12 show a comparison of the 
evolution of the two instabilities for \x = 0.5, v — 2, and (3 = —3. The system 
size is L = 104 and the numerical resolution is 256 x 256. The wavenumbers 
of the initial, slightly perturbed regular hexagon patterns are chosen care- 
fully to obtain the same linear growth rates for both modes. In particular, for 
q = —0.48 (q = 0.17) only the longitudinal (transverse) mode is destabilizing 
(cf. fig.4). 

Fig. 10 shows the early evolution of the instabilities in terms of the contour lines 
of the real and imaginary parts of the three amplitudes A t , i — 1,2, 3. While 
the longitudinal mode shown in the top three panels induces compressions and 
dilations of the pattern, the transverse mode shears the amplitudes. Both lead 
to the formation of defects. However, the defects induced by the longitudinal 
mode are typically aligned along the rolls and form where the bulges are 
maximal. 




Fig. 10. Zero contour lines of the amplitudes (A\ to ^3 from left to right) for (3 = —3 
and [i = 0.5. Solid lines are zero contours of the real parts, and the dashed lines are 
the imaginary parts. The top three panels are the longitudinal modes (q = —0.48) 
and the bottom three panels are the transverse modes (q = 0.17) at t = 100 in the 
simulations. Only a quarter of the domain is shown. 
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Fig. 11. Number of defects as a function of time for the longitudinal case (solid line) 
and transverse case (dashed line). 

The temporal evolution of the number of defects is shown in fig. 11. For both 
modes, we define t — when the first defect is formed. Both modes having 
the same growth rate it is not surprising that in both cases the defect number 
grows on a similar time scale. In fact, both reach roughly the same number of 
defects at about the same time t ~ 50. In both cases the subsequent ordering 
of the pattern appears to occur in two stages, characterized by an initial rapid 
annihilation of defect pairs and a later much slower phase. In the longitudinal 
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Fig. 12. Reconstructed temperature for the longitudinal mode (left panel, 
q = —0.48) and the transverse mode (right panel, q = 0.17) at t = 1750. 

case the defect number decreases in large steps, which are associated with the 
annihilation of a string of defects roughly aligned with the rolls, whereas no 
such steps are visible in the transverse case. The overall decay is substantially 
slower for the pattern induced by the transverse instability. This had also been 
found in the absence of mean flow [29] . 

Snapshots of the reconstructed patterns at the final time t = 1750 are shown in 
fig. 12. Only a quarter of the whole system is shown. While in the longitudinal 
the defect density case is very low and the defects are essentially isolated 
from each other, in the transverse case most of the defects are part of grain 
boundaries. 

To quantify the evolution of the amount of disorder in the orientation of the 
hexagons we determine the orientation of the local wavevector qj relative to 
the roll direction fij as defined by the angle af 



a, = arctan 



9i • T 3 



Qj 



W 



A, 



for j 



1, 2, 3. 



(35) 



In figure 13 we plot the probability distribution function (PDF) of the orien- 
tation angle a\ as a function of time for the cases shown in figures 10 and 11. 
Again, time t = is where the first dislocation occurs, and we have truncated 
the large peaks (around a\ = 0) at early times so that more structures can be 
discerned at later times. Similar PDFs are found for the other two amplitudes 
as well. For the longitudinal mode the PDF is centered around a,\ — from 
the beginning to the end, with the peak broadening around t = when the 
first few defects appear. Thus, for all times there is only a single domain and 
the hexagons remain essentially aligned with their initial orientation. In the 
transverse case, however, the initial peak at a\ = quickly decays and gives 
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Fig. 13. Evolution of the probability distribution of u\: the angle between k\ and 
h\ = x for both longitudinal (left) and transverse (right). 

way to two peaks of comparable size. This occurs around the time when the 
maximum number of defects is reached. The bi-modal PDF indicates that the 
transverse mode predominantly induces hexagons of two different orientations, 
which then co-exist for a long time. 

Hexagons vs. Rolls. 

In the absence of mean flow, the competition between uniform hexagons 
and rolls is governed by the energy difference between them. For a given q, 
hexagons are more stable than rolls for /i lower than a threshold value fith(o), 
at which the rolls and hexagons have the same free energy [29], while rolls are 
energetically favored above n t h- This boundary is indicated by the dashed line 
in figure 3. Below the dashed line, outside the stability balloon, the unsta- 
ble transverse and longitudinal modes grow and evolve towards a hexagon of 
different wavenumber. Above this line, rolls appear during the transients and 
eventually replace the hexagon. 

Due to the lack of a Lyapunov functional for eqs.(4,5), we resort to numerical 
simulations to locate the boundary between hexagon and roll in the presence 
of mean flow. Table 1 lists fi t h fo r q = ±0.6 for different values of (3. For 
example, for j3 = —0.1 (Prandtl number = 10) and q = 0.6 we find rolls as 
the final state at ji — 1.45 (figure 14(a)) while a mixed state of rolls and 
hexagons is found at /i = 1.425 (figure 14(b)). Thus, the transition value in /i 
at q — 0.6 for (3 = —0.1 is between 1.425 and 1.45. The enhanced instability 
of the transverse mode for q = 0.6 leads apparently to an earlier transition to 
rolls when the Prandtl number is decreased. At this point it is not clear why 
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Fig. 14. Reconstructed patterns from simulating equations (4,5) for (3 = —0.1 and 
q = 0.6. On the left, u = 1.425 at t = 2000, and u = 1.45 on the right panel at the 
same time. 
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Table 1 
Competition between rolls and hexagons for q = ±0.6 for different Prandtl numbers. 

the converse is not the case for q = —0.6. 
5.2 Effect of mean flow on motion of defects 



In this section we study the effect of the mean flow on the motion of a penta- 
hepta defect (PHD), which is a bound state of two dislocations in two of the 
three amplitudes. For large Prandtl number, where the mean flow is negligible, 
there exists a semi-closed form solution for stationary penta-hepta defects for 
hexagons at the band-center (q = 0) [33]. We study the general case with mean 
flow numerically by embedding two PHDs in the system as initial conditions 
and measure their velocity as a function of fi and q for fixed \i — 1 and v — 2. 

The numerical resolution is 256 x 256 in a system of size L = 400 and the time 
step is fixed at 0.1 for most of the results presented in this subsection. To satisfy 
the periodic boundary conditions, we place two PHDs of charges (0, +1, —1) 
and (0,-1, +1) in the computational domain, i.e. each PHD consists of two 
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Fig. 15. Defect velocity as a function of the coupling strength (3 for q = and 
H = 1. As in the case of rolls [49], for Prandtl number Pr > 0.5 the defect speed 
scales as ~ Pr -1 (cf. fig. 2). 

dislocations of opposite charges in the amplitudes A 2 and A 3 . We also apply 
a circular ramp at R — 0.4L, beyond which the phase is set to constant 
[32]. The interaction between pairs of PHD is characterized by the number 
N = X^ =1 5]5?, where 5- is the charge of the first (second) PHD in the 
jth amplitude [33]. In our simulations iV = —2 and the PHDs attract each 
other. Their interaction decreases with distance and we find that it becomes 
negligible for distances larger than 300 (for j3 = and q = the interaction- 
induced velocity is then below v = 0.001). Thus in the following, we place the 
two PHDs at least at a distance of 300 apart in the initial configuration for 
the velocity measurement of individual PHD. 

In the absence of mean flow, each independent PHD is found to move at a 
constant velocity, which vanishes at q = [32,33]. In the presence of mean 
flow, we also find that isolated defects move at a constant velocity. It is shown 
in figure 15 as a function of /3 for q = 0. The range of the linear scaling with 
respect to the coupling strength indicates that, for small /3, the contribution 
of the mean flow to the PHD velocity is purely additive via the advective 
term i/3A(f ■ V)Q and that the amplitudes A^ of the defect solution are only 
weakly affected by the mean flow. However, as \f3\ increases, the defect velocity 
deviates significantly from the linear scaling. This is analogous to the effect 
of mean flow on dislocations in roll pattern [50]. Also, for larger \f3\ the mean 
flow structure becomes distorted near the defect as shown in figure 16. The 
mean flow consists of two pairs of vortices, and is almost zero away from the 
defect. While for /3 = —0.2 (left panel) the two vortex pairs are of comparable 
strength, the pair on the right is much stronger than that on the left when (3 is 
decreased to /3 — —1.4 (right panel) This change occurs smoothly in f3. When 
the Prandtl number is decreased further so that /3 is below —2 the stability 
limit comes very close to the background wavenumber q = of the pattern 
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(b) 

Fig. 16. Mean flow structure around a PHD (q = 0): Panel (a) for /3 = —0.2 and 
panel(b) for (3 = —1.4. The solid lines (dash lines) are the zero contour lines for the 
real (imaginary) parts of the second amplitude are shown. 
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Fig. 17. Defect velocity as a function of wavenumber q for \i = 1, (3 = (dashed 
line) and j3 = —0.2 (solid line). 

and the PHD's triggers the formation of additional defects in their vicinity. 

The defect velocity also depends on the wavevectors <& of the three modes 
making up the hexagon pattern [32]. More specifically, within equations (4,5) 
it depends only on the projections qi-hi. Figure 17 shows the defect velocity as 
a function of q = qi -fii for (3 = (dashed line) and (3 = —0.2 (solid line). The 
mean flow shifts the defect velocity to more positive values for all q, implying 
that the wavenumber q st at which the defect remains motionless is shifted 
from q st = to negative values. For situations in which the evolution from 
disordered to more ordered patterns is dominated by defects this suggests that 
the wavenumber of the final state is in general not the critical wavenumber 
or that with the maximal growth rate but depends on the Prandtl number 
through the mean flow. 
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Fig. 18. Panel(a): Reconstructed hexagon pattern (q c = 4) at t ~ 14000 for 
(3 = —1.24 (Prandtl number ~ 1) and [i = 1.2. The initial condition was an or- 
dered, unstable hexagon pattern with q = 0.60. Panel(b): Mean flow of the state 
shown in (a). 

In various simulations of the nonlinear evolution of the instabilities of the 
hexagon pattern we found disordered states characterized by grain boundaries 
between domains of hexagons of different orientation that moved exceedingly 
slowly. An example of such a long-lived disordered state and the associated 
mean flow is shown in flg.l8a,b, where the spatial structure of such long-lived 
grain boundaries and the corresponding mean flow is depicted for /i = 1.2, 
q = 0.6 and (3 = —1.24. Here the resolution is 128 x 128 for a system size of 
L = 42. 

To identify the origin of these slow dynamics we performed simulations start- 
ing from patterns with two straight grain boundaries separating two hexagon 
patterns of different orientation rotated by tc/2 relative to each other. For 
certain magnitudes of the initial wavevectors the grain boundaries did not an- 
nihilate each other but persisted for an exceedingly long time. Fig. 19a shows 
such an initial condition with q? = 0.2. A contour plot of the histogram of the 
local wavevector of each of the three modes in the initial condition fig. 19a is 
shown in 20a. Solid lines pertain to q%, dashed lines to #2 , and dashed-dotted 
lines to q 3 . Despite the large magnitude of the reduced wavevector in the cen- 
ter domain the pattern is still linearly stable since the longitudinal component 
of the wavevector vanishes, g} -hj = 0. Within (4,5) only this projection enters 
the stability conditions. Either without (/3 = 0) and with (/3 = —2) mean flow 
the initial condition evolves to the patterns shown in figs. 19b and c, respec- 
tively, for t = 80, 000 and pi = 0.5. For (3 = simulation beyond t = 80, 000 
showed little difference in the spatial structure, and the defects seem to have 
reached asymptotic, immobile states at the end of the simulations. In the 
(3 = —2.0 case, although the spectra and the domain sizes remain more of 
less the same after t = 80, 000, defects exhibit persistent lateral motion along 
the grain boundaries (velocity ~4x 10 -4 ) throughout the simulation which 
continues beyond t = 160, 000 and lead to a slightly fluctuating shape of the 
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Fig. 19. Zero contour lines of the real and the imaginary part of A 2 - Panel (a) is 
the initial conditions, and panels (b) and (c) are snapshots at t = 80000 for (3 - 
and j3 = —2.0, respectively. 
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Fig. 20. Wave vector spectra of the complex amplitudes shown in figure 19. Panel 
(a) is for the initial conditions, and panels (b) and (c) are snapshots at t = 80000 
for (5 = and (5 = —2.0, respectively. 

domains. 

The histograms of the local wavevector of the final states depicted in fig.20(b,c) 
show that without mean flow all three wavevectors are essentially perpendic- 
ular to hi implying that individual defects would not move. In the simulation 
with mean flow the wavevectors are clearly not perpendicular to hi. Instead, 
q 2 ■ h 2 ~ —0.1 and q 3 ■ h 3 rs —0.1. In separate simulations we find that for 
f3 = —2 and \x = 0.5 this is the wavenumber at which individual PHD's do 
not move, q s t = —0.1. The histogram for the third wavevector qi, however, has 
peaks at (q x ,Q y ) = (—0.1,-0.05) and (—0.05,0.6). Thus, the two peaks have 
different projections onto hi and only one of them agrees with q st . This may be 
related to the fact that Ax has no dislocations in the grain boundary. In sep- 
arate simulations of individual PHD's with dislocations in A 2 and A 3 we find 
that the velocity of the PHD depends only weakly on the wavenumber of the 
defect-free component (Ai), and over the whole range —0.15 < q\ -hi < —0.05 
the PHDs are essentially motionless. Specifically, the magnitude of the velocity 
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is below 5 x 10 4 for \x = 0.5, (3 = —2, q\ ■ hi = —0.1, and q^ ■ h% = —0.1. 

With and without mean flow, the orientation of the initial and of the final 
wavevectors in the top and bottom domains are different implying that in 
these domains the pattern rotated until the projections of their wavevectors 
reached q st . This suggests that, more generally, in the ordering dynamics of 
hexagons starting from random initial conditions the orientation of a hexagon 
pattern within a given domain may rotate in a similar fashion and in the 
long-time dynamics the orientation of adjacent domains may predominantly 
be close to the stationarity condition for PHD's, i.e. the projections of their 
wavevectors onto a suitably chosen n have the same value and that value is 
close to q st . Of course, since our results are based on Ginzburg-Landau equa- 
tions they only apply to grain boundaries across which the orientation changes 
only by small amounts. Furthermore, in the truncation (4,5) the higher-order 
transverse derivatives V • fj have been neglected. They are expected to mod- 
ify the stationarity condition, since the defect velocity will also depend on the 
transverse component of the wavevector, <&-fj. The independence of the veloc- 
ity on the transverse component within (4,5) is related to the isotropy of the 
system. Thus, it is expected that at higher orders the condition q} • hi — q s t is 
replaced by a more complicated, but qualitatively similar condition suggesting 
similar behavior for grain boundaries. 



6 Conclusion 



The importance of mean flows has been noted in a wide range of pattern- 
forming systems. Various different types can be distinguished. In systems like 
binary- mixture convection [51] and in surface waves on liquids with small 
viscosity [52] they are driven by traveling waves. In other systems they corre- 
spond to a Goldstone mode as is the case in free-slip convection (e.g. [40,53]) or 
they arise from a conservation law, e.g. in systems with a deformable interface 
[54,55]. Then they constitute a separate dynamical variable and satisfy an evo- 
lution equation of their own. Depending on the symmetries of the system, the 
additional degree of freedom can introduce a host of new phenomena (e.g. [56- 
58]). In this paper we have studied the mean flow that is driven in Rayleigh- 
Benard convection by deformations of the convection pattern, which becomes 
relevant even close to onset for fluids with low Prandtl number [5,39,40]. In 
the realistic case of no-slip boundary conditions it does not constitute an addi- 
tional dynamical variable. Due to the incompressibility of the fluid it can only 
arise in three-dimensional systems, i.e. in two-dimensionally extended pat- 
terns, and introduces then a non-local interaction. Experimentally, its most 
striking signatures are the skew- varicose instability and the appearance of spi- 
ral defect chaos. The experimental observation of the skew- varicose instability 
and of transients of spiral-defect chaos in standing waves in vertically vibrated 
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granular media [22,23] suggests that a similar mean flow may also be relevant 
in that system. 

The focus of our weakly nonlinear analysis has been the impact of the mean 
flow on the stability and dynamics of hexagonal patterns. We have extended 
the usual three coupled Ginzburg-Landau equations for the description of 
hexagonal patterns by an equation for a vertical-vorticity mode in direct ex- 
tension of previous work on roll convection at small Prandtl numbers [39,40]. 

In general, there are two side-band instabilities that limit the band of sta- 
ble wavenumbers of hexagons. In the absence of mean flow it is always the 
longitudinal long-wave mode that is the relevant destabilizing mode imme- 
diately above the saddle-node bifurcation at which the hexagons come first 
into existence, while the transverse long-wave mode is the relevant mode for 
larger amplitudes. The Rayleigh number for the cross-over from one to the 
other mode depends on the cross-coupling coefficient v, but for realistic value 
it always occurs very close to the saddle-node bifurcation. The longitudinal 
mode is therefore poorly accessible in the absence of mean flow (e.g. [29]). We 
found that in the long-wave limit only one of the two phase modes, the trans- 
verse mode, couples to the mean flow and consequently only its stability limit 
depends on the Prandtl number. Specifically, for sufficiently small Prandtl 
numbers (~ 2.6 for v — 2) the transverse mode is stabilized for wavenumbers 
below the critical wavenumber to the extent that it is preempted by the longi- 
tudinal mode over the whole range of Rayleigh numbers from the saddle-node 
bifurcation to the transition to the mixed mode. In this regime studies of the 
difference between the two modes should be accessible in experiments since 
the type of instability encountered depends on whether the stability limits are 
crossed at the low- or the high-g side of the stability balloon. 

Our simulations of the nonlinear evolution of the instabilities indicate that, 
compared to the longitudinal instability, the transverse instability leads to a 
considerably larger number of penta-hepta defects and to more grain bound- 
aries separating patches of hexagons that are rotated with respect to each 
other. While indications of this were also found in the absence of the mean 
flow [29], the mean flow makes the distinction clear enough to make it worth 
addressing experimentally. To do so, hexagon patterns with a wavenumber 
away from the band center need to be initiated. Recently it has been shown 
that such initial conditions can, in fact, be prepared by a suitable localized 
heating of the fluid [37,50]. Given the striking difference in the transients aris- 
ing from the instabilities of the two modes it would be interesting to bring 
these techniques to bear in this system. 

Although the coupling to the mean flow makes the system non-variational no 
regime has been identified in which oscillatory instabilities are relevant. We 
also find that the stability limits are always determined by long-wave instabil- 
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ities and no additional instabilities at finite wavenumbers arise. Furthermore, 
to the order considered, the coupling to the mean flow vanishes at the band 
center. This implies that the pattern is always linearly stable there. 

As is the case in roll convection, the mean flow also affects the motion of de- 
fects. For the penta-hepta defects relevant in hexagonal patterns we find that 
similar to the case of rolls the wavenumber at which the defect is stationary 
is shifted to wavenumbers smaller than the critical one. For coarsening ex- 
periments starting from random initial conditions one may therefore expect 
that the eventual wavenumber of the ordered pattern may be reduced corre- 
spondingly. Our simulation suggest that the dependence of the defect velocity 
on the wave vector allows one to predict which grain boundaries have a par- 
ticularly long life time. Furthermore, a persistent drift of the defects in the 
grain boundary is also observed in the simulations. One may also expect that 
similar to the case of dislocations in roll patterns [38] the mean flow may allow 
two penta-hepta defects to form stable pairs if the background wavenumber is 
between the critical wavenumber and that corresponding to stationary defects. 
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A Derivation of the nonlinear phase equations 



Here we derive the nonlinear phase equations at the codimension-two point 
(q( ct \ R,( ct >) where both cr t and 07 are zero. Mainly, to obtain explicit expres- 
sions for (q( ct \ R^) we consider weak mean flow, (3 <C 1. (cf. eq.(30)). We 
rescale X = 5x, Y = 5y, and T = 5 A t. In this expansion 5 and (3 are two inde- 
pendent small parameters. Here we expand the amplitudes Aj = rje ignj ^ x ' y ^ + ' l ^ j 
as 



Tj = R { o t] + 5 2 r j2 + 5 4 r j4 + ■■■, ( A. 1 ) 

<t>j = $(<f>jo + 5 2 <j) j2 + 5 4 0i4 + •••), J = 1, 2, 3, (A.2) 

with Rq given in eq.(30). The mean flow amplitude Q is expanded accord- 
ingly in 5 2 

Q = 5 2 Q 2 + 5 4 Q 4 + ---. (A.3) 
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We substitute the above expansions into eqs.(4,5) and solve them at successive 
orders for 5. The mean flow feeds back to the equations for the phases <pj via 
the amplitude modulations (cf. eq.(5)). Thus, at each order we first solve for 
the amplitudes and the mean flow in terms of the phases that were determined 
at the previous orders, and substitute these solutions into the phase equations 
to obtain the phases at the next order. At 0(5 2 ), up to first order in j3, 



r i2 = — /_ , . ; dx<i>x ~ o../i , — i d x<t>x - 3d y 4> y )p, (A.4) 



^22 



1 n , 3_ 

2(1 + 1/) 8i/(l + i/) 

(d x - V3d y )((j) x - V3(j)y 



1 



4j2(l + i/) 



^—[(2d x + y/3d y )<f> x + y/3d x <j> y ]P, (A.5) 

(Ox + VSdy)((j) x + V3(j)y) 




+ Q n , A (2d x - Vsdy)^ - V3d x cf> y }P, (A.6) 

oi/( 1 + I/) 

ir ln ^ K-Vx4 (A.7) 

4za/2 1 + v) 16i/ i 1 + ^J 



At cubic order we recover eq.(24). In order to go on to fifth order we require 
that both 07 and a t vanish (up to second order in f3). As in section 3, the 
solutions are expressed in terms of the translation modes <p x = — (020 + 03o) 
and 4> y = (020 - 030 Vv^- 

Repeating the same procedures at 0(5 A ), we obtain r^. The expressions are 
too long to be displayed here. Also, at this order it is impossible to solve for 
Q4 in closed form. We therefore take the Laplacian of the phase equations at 
this order and substitute r,- 4 and V 2 Q4 to obtain the nonlinear equations for 
<f> x and (fi y at the codimension-two point: 

d t V 2 (f> x = £ x + NL x0 + pNL xl , (A.8) 

d t V 2 0y = A/ + NL ?,o + /3NL yl , ( A.9) 

where C x and C y are the linear terms given on the right-hand-sides of eq.(31) 
and eq.(32), respectively. The nonlinear terms at zeroth order in /3 for both 
equations, denoted as NL^o and NL^o, read as follows: 



,2 



NL, = -J^V 2 



is) 2 



d x <p x ( ^±^0, + ^^ + (3 + Su)d^ y ) , 
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3 + 5z/ o , 3(3 + v , ,. } 
dy<t>y ( " V " ' "' %<l>x + ^^T^d 2 ^ + 3(3 + 2v)d' xy < ,,„ 



{d y <p x + 8 x <p y ) ( (3 + 5i/)^ x + —2— ^ + -L-^^ 1 + 

/3 ( 3 + ^ , 3 + n ^2. , ,„ , 0,^2^ 



(A.10) 



NUn = -7^4v 



<?0^' ^ 2 



<w /3 + 8z/ o2 , 3 + 5i/_~ 9(1 + 1/) -a , \ 

d x <t> x -^—d%<t> x + — 2~ d 2 x <p Y + 2 ' ffiJ + 



{dy^ + d x <j> y ) —,—dl^ + -L-^tfU + 3(3 + ^ . + 



, 2 - , 3(3 + ^)^. , 27(l + z/) A2 



^ 3(3 + 2i/))O x + -^-^0, + v n y a : 



, (a.ii; 



with go defined in eq.(30). The expressions for NL^i and NL^i are too long to 
be shown here. 
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